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1. Project Summary 

Detailed quantitative spectroscopy of Type la supernovae (SNe la) provides crucial 
information needed to minimize systematic effects in both ongoing SNe la observational 
programs such as the Nearby Supernova Factory, ESSENCE, and the SuperNova Legacy 
Survey (SNLS) and in proposed JDEM missions such as SNAP, JEDI, and DESTINY. 

Quantitative spectroscopy is mandatory to quantify and understand the observational 
strategy of comparing "like versus like". It allows us to explore evolutionary effects, from 
variations in progenitor metallicity to variations in progenitor age, to variations in dust with 
cosmological epoch. It also allows us to interpret and quantify the effects of asphericity, as 
well as different amounts of mixing in the thermonuclear explosion. 

While all proposed cosmological measurements will be based on empirical calibrations, 
these calibrations must be interpreted and evaluated in terms of theoretical explosion models. 
Here quantitative spectroscopy is required, since explosion models can only be tested in detail 
by direct comparison of detailed NLTE synthetic spectra with observed spectra. 

Additionally, SNe IIP can be used as complementary cosmological probes via the spec- 
tral fitting expanding atmosphere method (SEAM) that we have developed. The SEAM 
method in principle can be used for distance determinations to much higher z than Type la 
supernovae. 
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We intend to model in detail the current, rapidly growing, database of SNe la and 
SNe IIP. Much of the data is immediately available in our public spectral and photometric 
database SUSPECT, which is widely used throughout the astronomical community. 

We bring to this effort a variety of complementary synthetic spectra modeling capabil- 
ities: the fast parameterized 1-D code SYNOW; BRUTE, a 3-D Monte-Carlo with similar 
assumptions to SYNOW; a 3-D Monte-Carlo spectropolarimetry code, SYNPOL; and the 
generalized full NLTE, fully relativistic stellar atmosphere code PHOENIX (which is being 
generalized to 3-D). 

2. Cosmology from Supernovae 

While indirect evidence for the cosmological acceleration can be deduced from a combi- 
nation of studies of the cosmic microwave background and large scale structure (Efstathiou 
et al. 2002; Bennett et al. 2003; Eisenstein et al. 2005), distance measurements to supernovae 
provide a valuable direct and model independent tracer of the evolution of the expansion 
scale factor necessary to constrain the nature of the proposed dark energy. The mystery 
of dark energy lies at the crossroads of astronomy and fundamental physics: the former is 
tasked with measuring its properties and the latter with explaining its origin. 

Presently, supernova measurements of the cosmological parameters are no longer limited 
by statistical uncertainties, but systematic uncertainties are the dominant source of error (see 
Knop et al. 2003, for a recent analysis). These include the effects of evolution (do SNe la 
behave in the same way in the early universe?) , the effect of intergalactic dust on the apparent 
brightness of the SNe la, and knowledge of the spectral energy distribution as a function of 
light curve phase (especially in the UV where are current data sets are quite limited). 

Recently major ground-based observational programs have begun: the Nearby Super- 
Nova Factory (see Aldering 2000; Nugent &: Aldering 2000), the European Supernova Re- 
search Training Network (RTN), the Carnegie Supernova Project (CSP), ESSENCE, and 
the SuperNova Legacy Survey. Their goals are to improve our understanding of the utility 
of Type la supernovae for cosmological measurements by refining the nearby Hubble dia- 
gram, and to make the first definitive measurement of the equation of state of the imiverse 
using z < 1 supernovae. Many new programs have recently been undertaken to probe the 
rest-frame UV region at moderate z, providing sensitivity to metallicity and imtested pro- 
genitor physics. SNLS has found striking diversity in the UV behavior that is not correlated 
with the normal light curve stretch parameter. As precise knowledge of the X-correction is 
needed to use SNe la to trace the deceleration expected beyond z —1 (Riess et al. 2004), 
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understanding the nature of this diversity is crucial in the quest for measuring dark energy. 
We plan to undertake an extensive theoretical program, which leverages our participation 
with both SNLS and the Supernova Factory, in order to refine our physical understanding 
of supernovae (both Type la and II) and the potential systematics involved in their use as 
cosmological probes for the Joint Dark Energy Mission (JDEM). 

In addition to SNe la, the Nearby Supernova Factory will observe scores of Type IIP 
supernovae in the Hubble Flow. These supernovae will provide us with a perfect laboratory 
to probe the mechanisms behind these core-collapse events, the energetics of the explosion, 
asymmetries in the explosion event and thereby provide us with an independent tool for 
precision measurements of the cosmological parameters. 

The SEAM method has shown that accurate distances may be obtained to SXe IIP, even 
when the classical expanding photosphere method fails (see Fig. 1 and Baron et al. 2004). 
Another part of the SN IIP study is based a correlation between the absolute brightness of 
SNe IIP and the expansion velocities derived from the Fe II 5169 A P-Cygni feature observed 
during their plateau phases (Hamuy &: Pinto 2002). We have refined this method in two 
ways (P. Nugent et al, 2005, in preparation) and have applied it to five SNe IIP at z < 0.35. 
Improving the accuracy of measuring distances to SNe IIP has potential benefits well beyond 
a systematically independent measurement of the cosmological parameters based on SNe la 
or other methods. Several plausible models for the time evolution of the dark energy require 
distance measures to 2; ~ 2 and beyond. At such large redshifts both weak lensing and SXe la 
may become ineffective probes, the latter due to the drop-off in rates suggested by recent 
work (Strolger et al. 2004). Current models for the cosmic star-formation history predict 
an abundant source of core-collapse at these epochs and future facilities, such as JDEM, in 
concert with the James Webb Space Telescope (JWST) or the Thirty Meter Telescope, could 
potentially use SNe IIP to determine distances at these epochs. 

Spectrum synthesis computations provide the only way to study this wealth of data and 
use it to quantify and correct for potential systematics and improve the distances measure- 
ments to both SNe la and SNe IIP. 



3. Understanding the 3-D Nature of Supernovae 

While most SNe la do not show signs of polarization, a subset of them do. These su- 
pernovae will play a role in determining the underlying progenitor systems/ explosion mech- 
anisms for SNe la which is key to ascertaining potential evolutionary effects with redshift. 
Flux and polarization measurements of SN 2001el (Wang et al. 2002) clearly showed polar- 
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ization across the high-velocity Ca II IR triplet. A 3-D spectopolometric model fit for this 
object assumes that there is a blob of calcium at high-velocity over an ellipsoidal atmosphere 
with an asphericity of ~ 15% (see Fig 3 and Kasen et al. 2003). Kasen & Plewa (2005) have 
shown that a gravitationally confined thermonuclear supernova model can also explain this 
polarization signature. If this is in fact the correct hydrodynamical explosion model for 
SNe la, then the parameter space for potential systematics becomes significantly smaller in 
their use as standard candles. Currently there are a wide variety of possible mechanisms to 
make a SN la each with its own set of potential evolutionary systematics. Thomas et al. 
(2004) showed that the observed spectral homogeneity implies that arbitrary asymmetries 
in SNe la are ruled out. The only way to test detailed hydrodynamical models of the ex- 
plosion event is to confront observations such as those that will be obtained via the Nearby 
Supernova Factory with the models via spectrum synthesis. The importance of studying these 
events in 3-D is clear from the observations, and therefore every effort must be made to 
achieve this goal. 



4. PHOENIX 

In order to model astrophysical plasmas under a variety of conditions, including differ- 
ential expansion at relativistic velocities found in supernovae. we have developed a powerful 
set of working computational tools which includes the fully relativistic, non-local thermody- 
namic equilibrium (NLTE) general stellar atmosphere and spectral synthesis code PHOENIX 
(Hauschildt k Baron 2004, 1999; Hauschildt et al. 1997a,b; Allard et al. 1997). PHOENIX is 
a state-of-the-art model atmosphere spectrum synthesis code which has been developed and 
maintained by some of us to tackle science problems ranging from the atmospheres of brown 
dwarfs, cool stars, novae and supernovae to active galactic nuclei and extra-solar planets. 
We solve the fully relativistic radiative transport equation for a variety of spatial boundary 
conditions in both spherical and plane-parallel geometries for both continuum and line radi- 
ation simultaneously and self-consistently. We also solve the full multi-level NLTE transfer 
and rate equations for a large number of atomic species, including non-thermal processes. 

To illustrate the nature that our future research will take, we now describe some of the 
past SN la work with PHOENIX. Nugent et al. (1995), showed that the diversity in the peak 
of the light curves of SNe la was correlated with the effective temperature and likely the 
nickel mass (see Fig. 2). We also showed that the spectroscopic features of Si II and Ca II 
near maximum light correlate with the peak brightness of the SN la and that the spectrum 
synthesis models by PHOENIX were nicely able to reproduce this effect. We were able to 
define two spectroscopic indices and (see Figs 4-5), which correlate very well with 
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the light curve shape parameter (Garnavich et al. 2004). These spectroscopic indices offer 
an independent (and since they are intrinsic, they are also reddening independent) approach 
to determining peak luminosities of SNe la. S. Bongard et al. (in preparation) have shown 
that measuring these spectroscopic indicators may be automated, and that they can be used 
with the spectral signal to noise and binning planned for the JDEM missions SNAP and 
JEDI. 

The relationship between the width (and hence risetime) of the lightcurves of SXe la to 
the brightness at maximum light is crucial for precision cosmology. It is well known that the 
square of the time difference between explosion and peak brightness, t^jg^ is proportional to 
the opacity, k, (Arnett 1982; Baron et al. 1993). In an effort to find a more direct connection 
between SN la models and the light-curve shape relationships we examined the Rosseland 
mean opacity, n, at the center of each model. We found that in our hotter, more luminous 
models k was a factor of 2 times greater than in our cooler, fainter models. This factor 
of 1.4 in trisc is very near to what one would expect, given the available photometric data, 
for the ratio of the light-curve shapes between the extremes of SN 1991T (an over-luminous 
SN la with a broad light curve) and SN 1991bg (an under-luminous SN la with a narrow 
light curve). 

We have been studying the effects of evolution on the spectra of SNe la, in particular 
the role the initial metallicity of the progenitor plays in the peak brightness of the SN la. 
Due to the effects of metal line blanketing one expects that the metallicity of the progenitor 
has a strong influence on the UV spectrum Hoflich et al. (1998); Lentz et al. (2000). In Lentz 
et al. (2000) we quantified these effects by varying the metallicity in the unburned layers and 
computing their resultant spectra at maximum light. 

Finally we note the work we have done on testing detailed hydrodynamical models of 
SNe la (Nugent et al. 1997). It is clear from these calculations that the sub-Chandrasekhar 
"helium-igniter" models (see for example Woosley & Weaver 1994) are too blue in general 
and that very few spectroscopic features match the observed spectrum. On the other hand, 
the Chandrasekhar-mass model W7 of Nomoto et al. (1984) is a fairly good match to the 
early spectra (which are most important for cosmological applications) of the most typical 
SNe la. Lentz et al. (2001) calculated an extensive time series of W7 and compared it with 
that of the well observed nearby SN la SX 1994D. In this work we showed that W7 fits 
the observations pretty well at early times, but the quality of the fits degrades by about 
15 days past maximum light. We speculate that this implies that the outer layers (seen 
earliest) of W7 reasonable well represent normal SNe la, whereas the inner layers of SNe la 
are affected by 3-D mixing effects. With the work described here, we will be able to directly 
test this hypothesis by calculating the spectra of full 3-D hydrodynamical calculations now 
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being performed by Gamezo and collaborators and by the Munich group (Hillebrandt and 
collaborators). Baron et al. (2005) have calculated very detailed NLTE models of W7 and 
delayed detonation models of Hoflich et al. (2002) . We find that W7 does not fit the observed 
Si II feature very well, although it does a good job in other parts of the spectrum. The 
delayed-detonation models do a bit better, but a highly parameterized model is the best. 
We will continue this work as well as extending it to 3-D models. This will significantly 
impact our understanding of SNe la progenitor, something that is crucial for the success of 
JDEM. 

We stress that the quantitative spectroscopic studies discussed here do not just show 
that a proposed explosion model fits or doesn't fit observed spectra, but provides important 
information into just how the spectrum forms. One learns as much from spectra that don't 
fit as from ones that do. 

Our theoretical work provides important constraints on the science definition of JDEM, 
helps to interpret the data coming in now from both nearby and mid-redshift surveys and 
involves ongoing code development to test 3-D hydrodynamic models, as well as both flux 
and polarization spectra from nearby supernovae which may indicate evidence of aspheric- 
ity. Research such as this requires both manpower and large-scale computational facilities 
for production which can be done to some extent at national facilities such as the National 
Energy Research Supercomputing Center at LBNL (NERSC), and local, mid-sized comput- 
ing facilities for code development which requires with the ability to perform tests with 
immediate turn-around. 
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Fig. 1. — The synthetic spectra of SN 1999em (dashed hnes) are compared to observed 
spectra (sohd hnes) at 5 different epochs. The abihty to accurately model the SED allows 
us to obtain accurate distances via the SEAM method (Baron et al. 2004). 
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Fig. 2. — The synthetic spectral sequence of SNe la compared with the observed sequence 
(Nugent et ah 1995). The existence of this spectral sequence and its correlation to the 
photometric one is an ideal example of what can be learned via quantitative spectroscopy 
that is much deeper than "this hydrodynamical model doesn't fit". 
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Fig. 3. — Results of 3D spectrum synthesis using SYNPOL for the SN la 2001el. Shown are 
the fits to the high- velocity Ca II IR triplet for various lines of sight assuming a clump of cal- 
cium obscures an elliptical photosphere (see Kasen et al. 2003; Kasen & Plewa 2005). Given 
enough observations at early times that show high-velocity calcium, we can test whether this 
model is correct and this leads to clues to both the progenitor system, explosion models, and 
potential systematics at high redshift. 
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Fig. 4. — The automated construction used for defining ^si- Green lines locate the Si ll A6355 
and Si ll A5972 Si ll lines. The blue lines locate the regions where the minima are searched 
for. The red diamonds locate the three maxima used to draw the two purple reference lines, 
and the two purple arrows represent d^ed and (ibiue used to compute = 
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Fig. 5. — The automated construction used to define ^ca- Red lines locate the 3650A and 
the 3933A Ca II lines. The maxima used to calculate are searched in the two blue 
colored regions. 



